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ABSTRACT 

Supermassive black holes are thought to be relics of quasars, and their numbers and 
masses are therefore related to the quasar luminosity function and its evolution with redshift. 
We have used the relationship between black hole mass and bulge velocity dispersion (the 
M. — a relation) to make an improved estimate of the mass density and mass spectrum of 
supermassive black holes. Uncertainties in the M. — a relation have little effect on the mass 
density. We find a mass density of (4.8 ± 1.6)/i 2 • lO 5 M Mpc" 3 . Some of the variance m 
published density estimates comes from the use of different values of the Hubble constant. 

Subject headings: galaxies: general, nuclei, black holes 



1. Introduction 

It is now fairly well established that most or all bulges contain supermassive black holes (Richstone 
et al. 1998; Kormendy & Gebhardt 2002; Magorrian 1998). Within the last two years it has become 
clear that the mass of the black hole (BH) is related to the velocity dispersion of the bulge of the host 
galaxy. This relationship has less scatter than a previously noted relationship between the black hole 
mass and the bulge light (Gebhardt et al. 2000; Ferrarese & Merritt 2000; Kormendy & Gebhardt 
2002). It has also been known for some time that the masses, and cosmic mass density of BHs, 
crudely correspond to the masses and numbers predicted by the luminous output of quasars, under 
the assumption that they are powered by accreting black holes with radiative efficiencies near 10%. 
Published estimates of the BH mass density predicted by quasars are around 2 • 10 5 M & Mpc~ 3 (Chokshi 
& Turner 1992; Soltan 1982; Salucci et al. 1999), while estimates from the X-ray background tend to 
be larger. Using the X-ray background, Salucci et al. (1999) obtain an estimate of 3 — 5 ■ l0 5 M & Mpc~ 3 , 
while Fabian & Iwasawa (1999) estimate the mass density to be 6 — 9 • 10 5 M Q Mpc~ 3 . 

It is useful to note that the above density estimates are independent of the Hubble constant or 
cosmological model (so long as it is general relativistic, isotropic and homogeneous) because they can 
be cast in terms of the surface brightness of quasar light on the sky, and hence determine its local 
energy density. Since the mass density of the radiating objects declines as the universe expands as 
(1 + z) 3 and the photon energy of their radiation declines as (1 + z) 4 , the only cosmological correction 
connecting the observed radiation to the local BH density is (1 + z), where z is the redshift of the 
epoch of the emission of the radiation. 

On the other hand, the local density of supermassive BHs are obtained by multiplying the density 
of galaxies by the typical BH mass per galaxy. Kormendy & Richstone (1995) used the BH mass-bulge 
luminosity relationship, together with the Loveday et al. (1992) estimate of the luminosity density 
in galaxies, to estimate the local density of BHs; for a Hubble constant of 80 km/s/Mpc this is 
1 • l0 6 M Q Mpc~ 3 . More recently, Merritt & Ferrarese (2001a) estimated this same quantity using the 
newer M — a relation. Their approach was to use the M — a relation to calibrate a mass-luminosity 
relation and then multiply by the galaxy luminosity function. 

In this paper we take an alternate approach. We use multiple estimates of the luminosity functions 
of galaxies for different Hubble types, together with estimates of bulge luminosity as a function of 
Hubble type, together with observed relationships between luminosity and velocity dispersions or 



- 3- 



rotation velocity, to construct a 'velocity function' (the number of galaxies per magnitude per unit 
volume) for bulges as a function of Hubble type. This function can then be used together with the 
M — a relation to estimate BH mass functions. Our basic result using this procedure is that the mean 
mass density of supermassive BHs is A.8h W0 2 ■ 10 5 M Q Mpc~ 3 . 

After we completed much of the analysis in this paper, we became aware of an estimate by 
Yu & Tremaine (2002) based on the Sloan Digital Sky Survey (SDSS), that gives an estimate of 
p, = (2.5 ± 0.4)/i .65 2 • W 5 M Q Mpc~ 3 . By assuming that h=0.65, as in their paper, our results are 
comparable with theirs. The rather surprising conclusions about the radiative efficiency of quasars 
derived by Yu & Tremaine (2002) from their low value of p. illustrate the importance of an improved 
estimate of this parameter. 



2. Method 

The basic assumption underlying this note is that the bulge velocity dispersion of a galaxy can 
be used to predict the mass of its supermassive black hole. Under that assumption, the best way to 
proceed would be to directly observe velocity dispersions of a volume-limited sample of galaxy bulges 
and to integrate over the implied distribution. Such a catalogue might well exist in the future, 3 but at 
present we adopt a simpler expediency. We use a power law (Faber & Jackson 1976) to convert various 
luminosity functions to a "dispersion function" , and the M — a relation for black holes to convert the 
dispersion function to a BH mass function. 

In general, if there are N galaxies per cubic megaparsec brighter than a luminosity L, then the 
luminosity function is 

This can be directly converted to the BH mass function by using the derivative of the relationship 
between luminosity and the velocity dispersion (a) and the black hole mass (A4») and a for bulges. 
Without any loss of generality, it can be written as 

dN ^ /rx fdL\ (dMS " J 



Since in what follows we will rely on luminosity functions due to Marzke et al. (1994), we adopt 
his parameterization of the luminosity function of different types of galaxies as Schechter functions and 
adopt power laws for the two relationships. To simplify the algebra we normalize luminosity, dispersion 
and black hole mass by their values at the fiducial luminosity of the Schechter function L*, as follows 

L = LjU (3) 
a = a/a* (4) 
M. =M./M„ (5) 



3 After starting this project, the early release data for the SDSS became available, but we pursued 
this calculation from completely independent data. 
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We then can write the Schechter parameterization and the relationships between luminosity and 
dispersion and dispersion and A4. as 

% = *W = (fe) ^ (6) 
L = a n (7) 
M. = a x (8) 

with 0* in Mpc~ 3 . Substituting these equations into equation 2 gives the following dispersion function 

dN_f &>-i e -*" ( 9 ) 



da 

with j3 = n(a + 1) and the following mass function for supermassive black holes 

-^L = e(M.) = e^M^M. 1 ' 1 e- M - e (io) 

dM. 

with = e0*, e = n/X and 7 = (5/X — n(a + 1)/A. 

Unless explicitly parametrized by h = H Q / (100km ■ s' 1 ■ Mpc' 1 ), all calculations were done 
assuming H Q = 100km ■ s" 1 • Mpc^ 1 . 



3. Results 

In order to determine the number density as a function of dispersion, the Schechter function fits, 
valid for M z > —21.5 (corresponding to a* = 400 km/s and M.,* = 1.8 • 10 9 M Q for an elliptical 
galaxy), obtained by Marzke et al. (1994), were modified using the corresponding bulge to total 
galaxy luminosity ratios, calculated from the data given in Simien & de Vaucouleurs (1986) for the 
appropriate galaxy range, as given in Table 1, to produce a fit as a function of the bulge-only luminosity 
for the spiral galaxies: 



dN 



dL 



bulge 



A power-law relation of the form 



*( W) = (fit) (ii) 

Lbulge ^bulgej ^- J *bulge (1^) 

L = a n (13) 



where L refers to the total luminosity ratio for the elliptical galaxies and to the bulge-only luminosity 
ratio for the spiral galaxies, may be used to convert the function of luminosity (6, 11) to one of 
dispersion such that 

^ - h 3 ■ (^) ~o^e-~° n (14) 



da 

with (5 = n(a + 1). 

To obtain a value for the parameter a* for each galaxy type, the parameter M* from the Marzke 
et al. (1994) Schechter function fits are converted to the equivalent dispersion, a* using a relationship 
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between the magnitude and the dispersion given by de Vaucouleurs & Olson (1982). The "best" 
relationship between magnitude and dispersion, selected from several in this paper, was determined by 
Gonzalez et al. (2000) to be: 

-{M Bt + hlog(h)) = a + b(log{a) - 2.3) (15) 

with a = 19.71 ± 0.08 and b = 7.7 ± 0.7. This results in a slope for the power-law relation (13) of 
n = 6/2.5 = 3.08 ± 0.28. However, this equation is in the M Bt magnitude system, and so must be 
converted into the Zwicky magnitude system in order to be combined with the Marzke et al. (1994) 
Schechter function fitting parameters. Despite concerns that there may not be a linear relationship 
between these two magnitude systems (Marzke, Huchra, & Geller 1994), an approximate conversion 
was adopted by combining the conversion given in Shanks et al. (1984) 

M b . =M Z - 0.45 (16) 

with that given in Gonzalez et al. (2000) 

M bj « M Bt + 0.06 (17) 

to obtain 

M Bt ^M z - 0.51 (18) 
which allows the magnitude-dispersion relation to be written as 

-M z + 0.51 = a + b(log(a) - 2.3) + 5log(h) (19) 

This relation may then be combined with the Simien & de Vaucouleurs (1986) magnitude correction 
for the bulge-to-total galaxy magnitude ratio, AM, so that only the M* bu i ge is used to calculate the 
dispersion for the spiral galaxies, such that 

-(M z - 0.51 + AM) = a + b(log(a) - 2.3) + 5log(h) (20) 

Hence, the velocity dispersion, a*, may be expressed as 

= 199.5 • 10 x (21) 

with 

x = (-M z * + 0.51 -a- 5log(h) - AM)r' (22) 

The number density as a function of dispersion (14) may then be combined with the relation 
between the dispersion and the black hole mass given by Tremaine et al. (2002), corrected from h=0.8 
(to h=l) for consistency, 

M, = hT x ■ 0.8 • 10 813±0 06 • " (23) 

(with A = 4.02 ± 0.32) to obtain the final equation for the number density as a function of black hole 

mass: 

= Q(M.) = h 3 ■ e^M^Mr 1 e' M ' £ (24) 

dM, 

again with 6* = e0*, e = n/X = 0.8 ± 0.1, 7 = (5/X and 

M.* = h- 1 ■ 0.8 • 10 8 - 13±0 - 06 • 10 A ' X (25) 
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This method of calculating rather than directly applying a relationship between A4, — Lbui ge , 

was used because Gebhardt et al. (2000) indicate that there is a tighter correlation between M, — a 
than between M., — Lb u i ge . 

Using expression (24) for the number density, an expression for the total mass density of black 
holes may be obtained by integrating the function 

P =fM. min MMM.)dM. (26) 
= h2 -SZ. m J*M:e~^dM. (27) 

from an assumed minimum black hole mass to an infinite black hole mass for each galaxy type, 
and then summing the results to obtain an estimate for the total mass density of black holes. The 
calculation was done adopting a minimum black hole mass limit of 1O 6 M , which results in a total 
mass density for all galaxy types of 

p. = (4.8±1.6)/i 2 -lO 5 M Mp C - 3 (28) 

The results of this calculation are given in Table 2; along with the number density function fitting 
parameters for each galaxy type. The uncertainties and average values given for each of these 
parameters were derived using a Monte Carlo analysis, with the largest source of uncertainty stemming 
from the bulge-to-total galaxy luminosity ratios given in Table 1. Figures 1-4 show histograms of the 
resulting densities from the Monte Carlo calculation, which illustrate that the median value for the 
density tends to be slightly lower than the average value, but still well within the stated uncertainties 
of the density, with a relatively symmetric distribution of densities. 

The calculation of the mass density was then redone using exactly the same method, but utilizing 
completely different galaxy luminosity functions obtained from Madgwick et al. (2002). In Madgwick 
et al. (2002) the galaxies were divided into subsets based on spectral type, 77, rather than by Hubble 
T-type. The parameter 77 is determined based on the absorption and emission line strengths in the 
galaxy spectrum. In order to apply the method developed in this paper, it was necessary to convert 
the spectral types into T-types so that the appropriate bulge-to-total luminosity conversions could be 
applied. This was done by utilizing the data given in Figure 4 of Madgwick et al. (2002) which showed 
the correlation between 77 and morphological type, and assuming a relatively even distribution of 77 over 
T-types. For the luminosity function designated as Type 1 (77 < —1.4), an assumption was made that 
galaxies were evenly divided between elliptical, SO, and Sa galaxies for the purposes of this calculation. 

The resulting mass density obtained using the Madgwick et al. (2002) luminosity functions is 

p. = (6.9 ± lA)h 2 ■ 10 5 M Q Mpc- 3 (29) 

with the results for the individual spectral types given in Table 3. This is slightly larger than the 
result obtained using the Marzke functions; however, it is strongly dependent on the assumptions 
made regarding the division into morphological type, particularly of the Madgwick et al. (2002) Type 
1 galaxies. If is is assumed that all of the Type 1 galaxies are elliptical, the resulting mass density 
becomes (10.5 ± 2.0)h 2 ■ 10 5 M Q Mpc~ 3 , if it is assumed that they are all SO galaxies then the mass 
density becomes (6.0 ± 2.0)h 2 ■ 10 5 M e Mpc~ 3 , and if it is assumed that they are all Sa galaxies, then 
the mass density is calculated to be (4.0 ± 1.8)h 2 ■ l0 5 M G Mpc~ 3 . Considering the uncertainties in using 
the Madgwick et al. (2002) spectral type classifications for the purpose of applying a bulge-to-disk 
luminosity correction, these results are consistent with those obtained using Marzke et al. (1994). 
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Also of importance is the fact that Madgwick et al. (2002) specifically take into account issues of 
Malmquist bias when determining their luminosity functions. Marzke et al. (1994) make no mention 
of this. Therefore, it appears that this possible omission will not have a significant impact on the final 
density of black holes. 



In Table 2, as expected, the ellipticals have the largest central black hole mass, Al,*, values, while 
later-type galaxies have increasingly smaller M.,* values, corresponding to their smaller relative bulge 
sizes. Plots of the number function (Figure 5) and the mass function (Figure 6) illustrate that the 
contribution from the SO galaxies dominates at low M,, while the ellipticals dominate for high M. 
values. Overall, the strongest contribution to the total mass density is from SO galaxies, with important 
contributions from the elliptical and Sa-Sb galaxies, but a small contribution from the Sc-Sd galaxies. 

The total cosmic mass density of black holes of (4.8 ± 1.6)h 2 ■ 10 5 M Q Mpc~ 3 obtained from this 
calculation is relatively independent of several of the adopted parameters. As illustrated by Figure 7, 
decreasing the lower limit of integration in the calculation of the mass density does not significantly 
change the resulting mass density estimate. This is because the black holes located in small-bulge 
galaxies are, themselves, small, and there are not enough small-bulge galaxies to make up in number for 
this lack in mass. This is further illustrated by the flat slope in the plot of the mass function (Figure 
6). Similarly, altering the value of A, the power law index, in the relationship between the black hole 
mass, A4., and the velocity dispersion, a (eqns 5 & 8), has a small effect on the final density. This is 
because most of the contribution is in the range of 10 7 — 1O 8 M - the same region in which the M.* for 
this relation is located, and in which the relationship is determined best. A perturbation in the power 
law index for masses near this mass will have little effect on the final answer. Changing the zero-point, 
however, would have a roughly linear effect on the resulting mass density, although for the currently 
accepted relationships, the differences in the zero-point are small. Redoing the calculation using the 
best-estimate Merritt & Ferrarese (2001b) relation stating that 



instead of the Tremaine et al. (2002) relation (23), produces a result of (4.3 ± 1.6)h 2 ■ 10 5 M o MpcT 3 , 
which agrees well with our previous estimate. This estimate for the mass density also agrees, 
within the limits of the uncertainty, with that given in the Merritt & Ferrarese (2001b) paper of 
p. « 3 ■ lO 5 M MpcT 3 . 

The parameters with the most significant impact on the resulting density, other than the M — a 
relationship from de Vaucouleurs & Olson (1982) and the Schechter functions from Marzke et al. 
(1994), are the assumptions about the morphological distributions of galaxies and the resulting 
bulge-to-disk luminosity corrections, adopted from Simien & de Vaucouleurs (1986). This is illustrated 
by the results obtained using the galaxy luminosity functions from Madgwick et al. (2002) instead 
of Marzke et al. (1994). By assigning a different morphological distribution to the Madgwick et 
al. (2002) Type 1 luminosity function galaxies, the resulting mass density varies by a factor of 
approximately 2.5. Given this uncertainty in the mass density resulting from the calculation using the 
Madgwick et al. (2002) data, we are more confident in our estimates derived from the Marzke et al. 
(1994) luminosity functions. 



4. 



Discussion 




(30) 
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The Schechter function fits may be less applicable in the low-black-hole mass regime, particularly 
for M, < 1O 6 M , and to a lesser degree for galaxies with M, < 1O 7 M . While there is information on 
black holes down to 1O 6 M , in our own galaxy and M32, the Schechter function fits may not correctly 
apply to elliptical galaxies, may include two types of galaxies within one classification, or may mix 
exponential and de Vaucouleurs-law galaxies together. We assume that while the Schechter function 
fits for elliptical galaxies may be dominated by dE galaxies at low luminosities, the same Ai, — o 
relationship is still applicable; in our own galaxy and a small number of other galaxies with exponential 
profiles and measured black hole masses, this relationship appears to hold true. 

It must be noted that the resulting mass density of black holes may be an underestimate due 
to a possible bias in the relationship between the black hole mass and the galaxy velocity dispersion 
(pointed out to us by Dr. Scott Tremaine). This bias relates to the fact that for any given velocity 
dispersion, there is a range of black hole masses, and this distribution may not be symmetrical about 
the ridgeline. We have investigated this possible bias by assuming a log-normal distribution of black 
hole masses, with a disperson of 0.15 in the log space, and recomputing the black hole mass density. 
This results in a mass density of p. = (5.0 ± 2.1)/i 2 • 10 5 M Q Mpc~ 3 , which is 4 % larger than our 
original estimate. 

It is often useful to have a simple approximation available for crude calculations. While the 
number functions we have derived are not of the Schechter form, we can find approximating functions 
of the form 

f/A ' C-M. ", M ' (31) 



dM. 

For the number function based on the Marzke et al. (1994) data, the roll off at high mass is too 
gradual to be well-approximated by any function of the form (31). We illustrate a fair fit in Figure 8. 
This fit, which works best near M., ~ 10 8 M Q , where the information on the numbers of black holes 
is the most extensive, has the parameters (c, .M,*, a) = (1.3 • lO~ lo M _1 Mpc~ 3 , 7.9 • 10 7 M Q ,0.95). 
We can find a better fit to the number function produced by the Madgwick et al. (2002) 
luminosity function. It is shown in Figure 9. The parameters for this fit are (c, M. 91f) a) 
= (3.2 • 10~ u M^Mpc- 3 , 1.3 • 1O 8 M , 1.25). 

The recent paper by Yu & Tremaine (2002) utilizes luminosity and velocity dispersion functions 
from early-type galaxies obtained as part of the Sloan Digital Sky Survey to perform a similar 
calculation. They obtain a result of (2.5 ± 0.4)/i 065 2 • 10 5 M Q Mpc~ 3 . If we take h=0.65 to compare 
our result to Yu & Tremaine (2002), we obtain p, = (2.0 ± 0.7) • 10 5 M Q Mpc~ 3 , which is in reasonable 
agreement with their result, particularly when the possible 4 % bias is taken into account. 

We acknowledge support for proposals 7388, 8591, 9106 and 9105, provided by NASA through 
grants from the Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-265555. 
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Fig. 1. — Histogram showing the 1000 computed mass densities for elliptical galaxies (in M & Mpc~ 3 ) 
obtained using a Monte Carlo analysis with a minimum black hole mass limit of 10 6 M Q . Both the 
average and median densities are 1.5 • lO 5 M Mpc~ 3 . 
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Fig. 2. — Histogram showing the 1000 computed mass densities for SO galaxies (in M Q Mpc~ 3 ) obtained 
using a Monte Carlo analysis with a minimum black hole mass limit of 10 6 M Q . The average density is 
2.1 • 10 5 M Q Mpc~ 3 , while the median density is 1.9 • 10 5 M Q Mpc" 3 . 
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Fig. 3. — Histogram showing the 1000 computed mass densities for Sa-Sb galaxies (in M Q Mpc~ 3 ) 
obtained using a Monte Carlo analysis with a minimum black hole mass limit of 1O 6 M . The average 
density is 1.1 • 10 5 M Q Mpc~ 3 , while the median density is 7.8 • lO 4 M Mpc~ 3 . 
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Fig. 4. — Histogram showing the 1000 computed mass densities for Sc-Sd galaxies(in M & Mpc 3 ) ob- 
tained using a Monte Carlo analysis with a minimum black hole mass limit of 1O 6 M . The average 
density is 1 • 10 4 M Q Mpc~ 3 , while the median density is 6 • l0 3 M Q Mpc~ 3 . 
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Fig. 5.— The log of dN/dM. in M^Mpc' 3 function of the black hole mass, M.» (in M Q )for each 
galaxy type. There is no information for M., < 1O 6 M ; the plot extends to 1O 5 M only to illustrate 
the extrapolation. 
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Fig. 6. — The log of Ai, dN/oLM. in Mpc~ 3 as a function of the black hole mass, A4. (in M Q ) for each 
galaxy type. There is no information for .M. < 1O 6 M ; the plot extends to 1O 5 M only to illustrate 
the extrapolation. 
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Fig. 7. — The log of the calculated black hole mass density (in M & Mpc 3 ) as a function of the minimum 
black hole mass (in M ) used in the integration of equation (27). 
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Log (Black Hole Mass) 

Fig. 8.— The fit to the dN/dM. (in M _1 Mpc~ 3 ) data derived from the original Marzke et al. (1994) 
luminosity functions, with Ai, in M & . The data is depicted by the solid line, while the dashed line is the 
function fit, as given in equation 31, with fitting parameters (c, Ai,*, a) = (1.3 • 10~ 10 Mq" 1 Mpc~ 3 , 7.9 • 
1O 7 M ,O.95). 



-18- 




Log (Black Hole Mass) 

Fig. 9.— The fit to the dN/dM. (in M^Mpc' 3 ) data derived from the original Madgwick et al. (2002) 
luminosity functions with A4, in M Q . The data is depicted by the solid line, while the dashed line is the 
function fit, as given in equation 31, with fitting parameters (c, .M,*, a) = (3.2 ■ 10 _11 M Q _1 Mpc~ 3 , 1.3 • 
1O 8 M ,1.25). 
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Table 1. Fitting Parameters for Galaxies. 



Galaxy Type 




0* 


a 


T Types 


AM 






(l(T 3 Mpc- 3 ) 










Elliptical 


-19.23 


1.5 ± 0.4 


-0.85 


N/A 


0. 


SO 


-18.74 


7.6 ± 2.0 


-0.94 


-3 - 





0.64 ±0.30 


Sa-Sb 


-18.72 


8.7 ± 2.2 


-0.58 


1 - 


4 


1.46 ±0.56 


Sc-Sd 


-18.81 


4.4 ± 1.1 


-0.96 


5 - 


7 


2.86 ±0.59 
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Table 2. Computed Mass Density for A4, min = 1O 6 M using Marzke et al. (1994) Parameters. 



Galaxy Type 


(1O 7 M ) 


(10~ 3 Mpc- 3 ) 


7 


P 

(lO 5 M Mpc~ 3 ) 


Elliptical 


11.2 ±1.9 


1.2 ±0.3 


0.12 ±0.01 


1.5 ±0.5 


SO 


3.1 ± 1.3 


5.9 ±1.6 


0.046 ± 0.006 


2.1 ± 1.1 


Sa-Sb 


1.3 ±1.0 


6.7 ±1.8 


0.32 ±0.04 


1.1 ± 1.0 


Sc-Sd 


0.29 ±0.26 


3.4 ±0.9 


0.031 ±0.004 


0.1 ±0.1 



Note. - For all galaxy types e = n/X — 0.8 ± 0.1. 
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Table 3. Computed Mass Density for M, m in = 1O 6 M using Madgwick et al. (2002) Parameters. 



Madgwick Galaxy Type 


Assumed T-Type(s) 


AM 




P 








(1O 7 M ) 


(lO 5 M Mpcr 3 ) 


Type 1 (rj < -1.4) 


Elliptical 


0. 


10.0 ±1.8 


3.0 ±0.6 


Type 1 (r) < -1.4) 


-3 - 


0.64 ±0.30 


4.9 ±2.0 


1.5 ±0.6 


Type 1 (77 < -1.4) 


1 - 2 


1.19 ±0.51 


2.7 ± 1.7 


0.81 ±0.50 


Type 2 (-1.4 < < 1.1) 


1 - 5 


1.58 ±0.60 


1.8 ±1.3 


1.1 ±0.9 


Type 3 (1.1 < r] < 3.5) 


3 - 5 


1.85 ±0.59 


0.86 ±0.70 


0.37 ±0.36 


Type 4 (77 > 3.5) 


5 


2.47 ±0.45 


0.37 ±0.26 


0.07 ±0.07 



